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Abstract 

Using combinatorial arguments, we determine an upper bound on achievable rates of 
stabilizer codes used over the quantum erasure channel. This allows us to recover the no- 
cloning bound on the capacity of the quantum erasure channel, R < 1 — 2p, for stabilizer 
codes: we also derive an improved upper bound of the form R < 1 — 2p — D(p) with a 
function D(p) that stays positive for < p < 1/2 and for any family of stabilizer codes 
whose generators have weights bounded from above by a constant - low density stabilizer 
codes. 

We obtain an application to percolation theory for a family of self-dual tilings of the 
hyperbolic plane. We associate a family of low density stabilizer codes with appropriate 
finite quotients of these tilings. We then relate the probability of percolation to the 
probability of a decoding error for these codes on the quantum erasure channel. The 
application of our upper bound on achievable rates of low density stabilizer codes gives 
rise to an upper bound on the critical probability for these tilings. 

1 Introduction 

Low Density Parity Check (LDPC) codes are classical error-correcting codes originally intro- 
duced by Gallager |19| . They come with highly efficient local iterative decoding schemes, have 
been extensively studied and have proved very successful on a number of channels. Therefore, 
in the field of quantum communication and quantum computation, it is natural to look into 
the quantum analog of classical LDPC codes, which arguably are stabilizer codes with gen- 
erators of bounded weight. We will call such codes low density stabilizer codes or, following 
others, refer to them somewhat loosely as quantum LDPC codes. These include a number 
of constructions of locally decodable quantum codes with a topological connection, starting 
with Kitaev's celebrated toric code [26J, and other families among which surfaces codes [D], 
|37| . color codes [TJ, [8j, and other variants |35| . |14j . Various generalizations to the quantum 
setting of classical LDPC codes have also been proposed, e.g. |T], [2], [13J, [23J, |32j . 

In the present work we are interested in the performance of quantum LDPC codes over 
the quantum erasure channel. Our motivation is inspired by the classical setting, in which the 
systematic study of LDPC codes for the classical erasure channel has led to a better under- 
standing of the behaviour of LDPC codes for more complicated channels such as the binary 
symmetric channel and the Gaussian channel. This approach has hardly been attempted in 



the quantum setting and it is not unreasonable to hope for similar returns in the long run. The 
quantum erasure channel, besides being simpler than the more universal depolarizing channel, 
is also a realistic channel [21]. Its capacity is known, it is 1 — 2p pJJ, however almost nothing 
is known about the performance of quantum LDPC codes over this channel. We would like to 
gain some understanding as to what are the code characteristics needed to achieve capacity. 

We shall derive a bound on the achievable rate of stabilizer codes as a function of an upper 
bound on the weight of the generators of the stabilizer group. Equivalently, we will derive an 
upper bound on the decoding threshold on the erasure channel for quantum LDPC codes. This 
bound will yield the following result: any family of stabilizer codes that have stabilizer groups 
with generators of weight bounded by a constant, cannot achieve the capacity of the quantum 
erasure channel. This phenomenon is somewhat analogous to the classical setting |19j |10| 
where it is known that capacity achieving LDPC codes must have parity-check matrices with 
growing row weights. 

Our result has an unexpected application to percolation theory. Given an infinite edge- 
transitive graph, a random subgraph, called the open subgraph is considered where every edge 
is declared open, independently of the others, with probability p. The central question in 
percolation theory is the determination of the critical probability p c , which is the minimum 
value of p such that the open connected component of any given edge e is infinite with non 
zero probability. For most graphs, computing the critical probability exactly is usually quite 
difficult. 

We are interested in the critical probability of graphs that make up regular tilings of the 
hyperbolic plane. The connection with quantum erasure correcting codes is that quotients 
of the infinite tiling of the hyperbolic plane yield finite graphs (more precisely combinatorial 
surfaces) that define quantum LDPC codes (surface codes): selecting a random erasure pattern 
for the finite code is the same as selecting a random subgraph of the finite graph, and the 
non-correctable erasure event is very close to the percolation event on the infinite graph. This 
connection was observed for the toric code |17| . and this is why the erasure threshold of the 
toric code coincides with the critical probability p c = 1/2 in the square lattice. We shall derive 
a bound on the erasure decoding threshold for surface codes that will lead to an upper bound 
on the critical probability p c for hyperbolic tilings. To the best of our knowledge this is the 
sharpest presently known such upper bound. 



The paper is organized as follows. After introductory and background material presented in 
Section [2] we derive an upper bound on achievable rates of stabilizer codes in Section [3] There 
are two main results in this section. The first is Theorem 13.51 which states that achievable rates 
of stabilizer codes satisfy a bound of the form R < 1 — 2p — D(p) for a non-negative function 
D(p). This enables one to recover the capacity bound R < 1 — 2p for the quantum erasure 
channel for the particular case of stabilizer codes. It also enables one to derive improved bounds 
on achievable rates for particular classes of stabilizer codes. We derive such an improved 
bound in Theorem |3.8| for stabilizer codes with generators of bounded weight. In Section [4] 
we refine the above bound for a particular class of LDPC CSS codes that make up a family 
of surface codes. In Section [6] we apply the refined bound on achievable rates to percolation 
on hyperbolic lattices: the main result is Theorem 6.9 which is an upper bound on critical 
probabilities. Finally, an appendix regroups some technicalities necessary to complete formal 
proofs (Appendix |A[) . 
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2 Background 



A stabilizer code of parameters [[n, k]] is a subspace of dimension 2 k of the space = 
(C 2 )® n . It is defined as the set of fixed points of an abelian group of Pauli operators. Below 
we go over notation and definitions. This material is quite well-known but we have felt the 
need to highlight the properties that we need, in particular because we shall make extensive 
use of linear algebra. For a more precise description of quantum information and quantum 
error-correcting codes, see Nielsen and Chuang [28J with a binary point of view close to the 
one that we adopt here, or the article of Calderbank, Rains, Shor and Sloane |12| . for an F4 
point of view. 

2.1 Pauli groups 

A quantum bit or qubit is a vector of % = C 2 . It is the basic unit of quantum information. 
A sequence of n qubits lives in the space %® n . The classical Pauli operators form a basis of 
the space of operators on %® n . 

Denote by / the identity matrix of size 2, X = (±q), Z = (0-1) an d Y = iXZ. These 
operators satisfy the following relations: 

' x 2 = Y 2 = Z 2 = I, 
XY = -YX = iZ, 
' YZ = -ZY = -iX, 
ZX = -XZ = iY. 

The Pauli group V\ for one qubit is the group generated by the matrices: 

V\ = {±1, ±U, ±X, ±iX, ±Y, ±iY, ±Z, ±iZ} 

Remark that two different non-identity Pauli matrices always anti-commute. 

The Pauli group V n on n qubits is the multiplicative group of n-fold tensor products of 
errors of V\: 

V n = {^Ex ®E 2 ®---®E n I o = 0,l,2 or 3 and E { = I,X,Y or Z} 

The complex number i a is the phase of the Pauli operator. An important consequence of this 
construction is the fact that two Pauli errors either commute or anti-commute. More precisely, 
given two errors E and E' of V n , we have: 

EE' = {-l)H E ' E '^E'E, 

where f(E,E') is the number of components j such that Ej and E'- are two different non- 
identity Pauli matrices. For example, the operators I®X®Z and X®Y®Z in V3 anti-commute 
because they anti-commute only in the second position. This fact is at the origin of syndrome 
measurement. 

2.2 Stabilizer codes 

A stabilizer group S is a commutative subgroup of V n which doesn't contain —I. A stabilizer 
group is generated by a family of commuting Pauli operators: S =< S\, S%, . . . , S r >. Without 
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loss of generality, we can assume that these operators have phase 1. This ensure us that —I 
is not in S. 

Given S a stabilizer group of V n , the corresponding stabilizer code C(S) is defined as the 
set of fixed points of the subgroup S in H® 71 . Using the physical ket notation for vectors, we 
have: 

C(S) = meH® n \sm = myseS}. 

This subspace is not trivial by construction of a stabilizer group. The integer n is the length 
of the quantum code. Assume that S is generated by r generators: S =< Si, S2, ■ ■ ■ , S r > 
with Si G V n - We call the stabilizer matrix of C(S) the matrix H G M r ,n({I, X, Y, Z}) 
with the i-th row representing the generator Si. The coefficient Hjj is the j'-th component 
of Si. For example, the quantum code associated with the 3 commuting generators Si = 
(I ® X <g> Z <g) Y ® Z), S 2 = (Z <g> Z <g> X ® I ® Z), and S3 = (J ® Y <g) Y ® Y ® Z) is described 
by the following stabilizer matrix: 

(X Z I I z\ 
H= \z X X Y l\ 

\Y Y X Y Zj 

A stabilizer code is completely defined by its stabilizer matrix, though different stabilizer 
matrices can define the same group and therefore the same code. 

2.3 Syndrome of an error 

Given S =< S±, S2, ■ ■ ■ , S r > a stabilizer group of V n , assume that \ip) G C{S) is subjected to 
a Pauli error E G V n . The vector is corrupted to E\ip). To recover the original quantum 
state, we measure the syndrome to obtain information on the error. The syndrome of E G V n 
is a(E) = (<ti, (72, . . . , u r ) G F2 defined by: 

I if E and Si commute 
o~i = < 

I 1 if E 1 and Si anti-commute 

Given the corrupted quantum state E\ip) the syndrome of the error E can be measured. It 
satisfies a (EE') = a(E) + a(E'). The syndrome of an error s G S which has no effect on the 
quantum code is a(s) = 0. 

2.4 Minimum distance of a stabilizer code 

The phase i a of a Pauli error E G Vn does not play a role because we want to protect quantum 
states of C(S), that is vectors of C(S) defined up to multiplication by a non-zero complex 
number. Therefore, we will consider errors E G V n defined up to phases. In what follows, 
unless otherwise stated, the Pauli group will be the abelian quotient group: 

V n = V n /{±l,±i}. 

Given E,E' in the group V n , we will say that they commute if they commute in the original 
group V n - This misuse of language is not problematic because commutation doesn't depend 
on the phase. 

If we receive a quantum state E\ip), where \t/j) is in the quantum code C(S), we measure 
its syndrome a(E). We then apply to E\ip) an error E such that a(E) = a(E). After this 
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process the quantum state is EE\ip). It is corrupted by an error EE of syndrome 0, because 
a (EE) = cr(E) + cr(E) = 0. There are two types of error of zero syndrome. If EE is 
in S, it fixes the quantum code and we have recovered the quantum state. Otherwise, the 
original quantum state is probably lost. Errors of zero syndrome that are not in S are called 
undetectable or problematic errors. 

The above observation leads to the definition of the minimum distance d: 

d = mm{\E\ | E G P n \S, a(E) = 0}, 

where \E\ is the weight of E, it is the number of non-identity components of E. In other 
words, d is the minimum weight of a problematic error. The set of errors of syndrome in V n 
is frequently denoted by N(S), because it is the normalizer and the centralizer of the subgroup 
S in V n . Thus, the minimum distance is the minimum weight of an error of N(S)\S. 

2.5 Degeneracy 

An essential feature of quantum coding theory that sets it appart from classical coding is 
degeneracy. It allows for the same decoding procedure to correct a large number of different 
errors. More precisely, all the errors of a coset E.S can be corrected by the same error E. 
Indeed, assume that a state of the quantum code is corrupted by an error Es, where s € S. 
Then, after application of E, we recover the original quantum state because it is fixed by S. 
We have: EEs\ip} = To correct an error E G V n , it is sufficient to determine its coset 
E.S. 

2.6 The F 2 -vector space structure, rank and dimension 

The dimension of the quantum code C(S) is 2 fc , where k = n — rankS". It is the number of 
encoded qubits. The quantity rankS" is simply the rank of the abelian group S, i.e. the size 
of a minimum generating family. The rate of the quantum code is defined as R = k/n. 

The Pauli group V n can be seen as an F2-vector space of dimension 2n, by the isomorphism: 

9 : V n — »• ¥2 x = F^" 
X t .— > ( ei \0) 
Zi .— )• (0\e t ) 

where X{ is X on the i-th component and the identity on the other components. Errors Y{ 
and Zi are defined similarly. The image of Y{ = X; L Zi is Q(XjYi) = (ej|ej). For example, 
the operator / % X % Y % Z corresponds to the vector (0110|0011). For this F2-vector space 
structure, the addition of vectors in F^™ corresponds to componentwise multiplication of Pauli 
errors. 

By the isomorphism 9, subgroups of V n are sent onto F2-linear subspaces of F| n . The rank 
of a subgroup of V n is therefore also the dimension of the corresponding subspace. If H is a 
stabilizer matrix we will also write rankH to denote the rank of its row-space, equivalently 
the rank of the associated stabilizer group. Note that we may choose a stabilizer matrix with 
a larger number r of rows than its rank. 

We will find it convenient to keep the notation /, X, Y and Z for stabilizer matrices, but 
we stress the binary vector space structure that we will rely upon heavily in the next section. 
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With this vector space interpretation, the syndrome application: 



a : V n .— > F r 2 
2? — ► cr(£) 

can be regarded as an F2-linear map. 
2.7 The CSS construction 

One of the most popular ways of constructing quantum codes is the Calderbank, Shor and 
Steane (CSS) construction [TTJ [M]. A CSS code is a stabilizer code constructed from a 
stabilizer group S such that: 

S =< S\, S2, ■ ■ ■ , S rx , SV X +1> • • • j Srx+rz ^ 

where Si, S2, • • . , S rx are included in {/, X}® n and S rx +i, S rx +2, ■ ■ ■ , S rx +r z belong to {J, Z}® n . 
This simplifies commutation relations because two errors of {I, X}®" automatically commute 
and it is similar in {I, Z}® n . In this case the stabilizer matrix H is decomposed into two 
stabilizer matrices and H^. The matrix is composed of rx rows representing the 
stabilizers with coefficients in {/, X} and the matrix is composed of rz rows which define 
the stabilizers with coefficients in {/, Z}. 

Remark that the subgroup {/, X} is isomorphic to F2, thus we can write the matrix Hx 
as a binary matrix. The same remark is also valid for Hz- By this last isomorphism, rows of 
the matrices can be seen as binary vectors of length n and the commutation relation between 
a row of Hjf and a row Hz corresponds to the orthogonality of these binary rows in FJj. 

Finally, a CSS code can be defined from two binary matrices Hx G Mr x ,n(^2) an d 
Hz 6 ■Mr z ,n(^2) with orthogonality between rows of Hj and rows of Hz in F^. The number 
of thus encoded qubits is: 

k = n — rank — rank , 

because ranks of the binary matrices Hx and Hz coincide with ranks of the corresponding 
groups. Denote by Cx the classical code KerHx = {c G FJj, Hx *c = 0} and denote by Cz 
the code KerH^- The minimum distance of the quantum code is: 

d = mf{w{x) I x e C x \Cz- U C z \Ci}, 

where is the Hamming weight of a binary vector. 



Problematic errors. By the isomorphism of section 2.6, the error vector E can be seen as 
two simultaneous binary vectors, 0(E) = (Ex,Ez)- The error E has zero syndrome if and 
only if 

E x e C z and E z G C x (1) 
The error E is problematic if ([T]) holds together with the condition 

E x ?C X or E Z $C%. (2) 
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3 Capacity of the quantum erasure channel 



The capacity of a quantum channel is the highest rate of a family of quantum codes with an 
asymptotic zero error probability after decoding. Such a rate is called achievable. For the 
quantum erasure channel of erasure probability p, the capacity Q is 1 — 2p when p < 1/2 and 
it is zero above 1/2. The upper bound 

Q < 1 - 2p (3) 

comes from the no-cloning theorem, see for example [5]. Therefore, it doesn't rely on the 
quantum code structure. Since our purpose is to obtain improved capacity bounds for partic- 
ular families of codes, namely quantum LDPC codes, we need to derive capacity from the code 
structure: our first step is to express achievable rates of stabilizer codes over the quantum 
erasure channel, as a function of their stabilizer matrices. 

3.1 The quantum erasure channel 

The quantum erasure channel admits several equivalent definitions. See for example |21| . |20] . 
|30| . As a completely positive trace preserving map, it is given by: 

MM'— K1-pMM+p|2><2| 

where is a quantum state in T~L and the final state lives in C 3 = % ® C|2). The vector 
|2) is orthogonal to the space T~L, it corresponds to a lost qubit. In this paper, we will use the 
definition based on the Pauli operators which is well adapted to the stabilizer formalism. When 
we use the quantum erasure channel, each qubit is erased independently with probability p. 
An erased qubit is subjected to a random Pauli error I, X, Y or Z with equal probability 1/4 
and we know that this qubit is erased. 

This description of the quantum erasure channel can be deduced from the definition as a 
completely positive trace preserving map. Indeed, the orthogonality between |2) and T~L allows 
us to measure the erased qubit. After, we replace the lost qubit |2)(2| by a totally random 
qubit of density matrix 1/2. This random state is the original qubit subjected to a random 
error /, X, Y or Z with equal probability. Therefore, we recover the second definition. 

On n qubits, we denote by £ G the characteristic vector of the erased positions. Each 
component of the vector £ follows a Bernoulli distribution of probability p. That is, the 
probability of a given vector £ G F2 is p' £ '(l — p) n_ ' £ L The qubit in position i is lost if and 
only if £j = 1. In this case, the quantum state is subjected to a random Pauli error E G V ni 
which act trivially on the non-erased qubits: E{ = I if £j = 0. We write this condition Bc£ 
and will say that erasure £ covers the error E. Keep in mind that E is a Pauli operator with 
coefficients in {/, X, Y, Z} and £ is a binary vector, the shorthand notation E C £ expresses 
just that the support of E is included in the set of erased positions. Note finally that given 
£ 6 Frj, all errors E C £ occur with the same probability. 

An encoded quantum state is corrupted to a state E\ip) by a random error E for which 
we have the additional knowledge E C £. To recover the original quantum state, we compute 
the syndrome a G F£ and must deduce from the couple (£,<r) an error E C £. To correct 
the effect of E we apply E and the final state is EE\t/j). If the errors E and E are in the 
same coset modulo S, then EE is a stabilizer of the quantum code. Thus the final quantum 
state is the original state. When E is not equivalent to E, we will not, in general, recover the 
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quantum state. Note that in this case EE is a problematic error and EE C £. When this 
happens, i.e. when the erasure vector covers a problematic error, we will say that we have a 
non- correctable erasure: otherwise the erasure is correctable. 

Non-correctable erasures in the CSS case. From the characterization ([TJ and ^ of 
problematic errors, we obtain the simple characterisation of a non-correctable erasure in the 
CSS case. 

Proposition 3.1. Let H = (hz) ^ e ^ e stabilizer matrix of a CSS code, and let Cx and 

Cz be the corresponding classical binary codes. The erasure vector £ G F% is non-correctable 
if and only if there exists a binary vector v whose support is included in the support of £ and 
such that 

veC x \C^ or veC z \C x . 



3.2 An example of a non-correctable erasure 

As an example of the general case, consider the stabilizer code defined by the matrix: 

(1 X Z Y Z^ 
H= Z Z X I Z 
\I Y Y Y Zj 

If the erasure is £ = (0, 1, 1, 0, 0), there are 2 2 ' £ ' = 2 4 possible errors: 

{E G V n I E C £} = {X^Z^X^Z^ | a u b t G F 2 }, 

where the error Xi is the error with z-th component X and which is the identity outside i. 
The operator Z{ is defined similarly and Yi is the error X{Z{. 
Let us focus our attention on the "erased matrix" 



which is the submatrix of H whose columns are the columns indexed by the erased positions. 
It is natural to introduce this matrix because the syndrome of an error included in the erasure 
£ depends only on these columns. We remark that the third row of this matrix is the product 
of the first two rows. Thus, the syndrome u G F| of an error E C £ satisfies 113 = u\ + U2- It 
depends only on the first two rows of Hg. Therefore, there are 2 2 different syndrome values 
for errors E that are covered by the erasure. 

Now let us look at the remaining columns. The non-erased submatrix Hg is: 



Assume that E and E' are two errors included in £, which are in the same degeneracy class. 
That is, they differ by right-multiplication by an error s G S. The restriction of the error 
s = EE' to £ is the identity. The rank of this submatrix is rankHg = 2 because the first 
row and the third row are identical. Therefore, we have two possibilities for s: either s = J® 5 
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or s = Si S3 = There are two errors in each degeneracy class and four 

possible syndrome values: therefore, if there were no problematic error included in £ the total 
number of errors included in £ would be 2 x 4 = 2 3 , but we have seen that it actually is 2 4 . 
This erasure is not correctable. 



3.3 Two enumeration lemmas 

We will pursue the preceding approach. Our strategy is to determine the cardinalities of two 
sets of Pauli errors: 

• N(S) E = {EgN(S) I EC £}, 

recall that N(S) is the set of Pauli errors of syndrome 0. 

• S £ = {s 6 S I s C £}. 

We will use the submatrices introduced in the above example. 

The random submatrix Hg: Let H be a matrix of a stabilizer code. With an erasure 
£ G Frj, we associate the submatrix Hg of the stabilizer matrix H G ■Mr.ni'Pi) composed of 
the columns of the erased qubits. This is the submatrix of the columns of index i such that 
£, = 1. Similarly Hg is the matrix of the non-erased qubits, corresponding to the conjugate 
£ of £ defined by: £j = £j + 1. 

Lemma 3.2. Let S be a stabilizer group of matrix H G Ai rn . The set iV(S)g is an ¥2-vector 
space of dimension 2|£| — rankHg. 



Proof. The F2-linear structure of the Pauli group V n has been detailed in Section 2.6 The 
syndrome is an F2-linear map from V n to F£. Its restriction erg to the space of Pauli errors 
included in £ is also an F2-linear map. The subspace iV(S)g is simply the kernel of erg. Its 
dimension is 2|£| — dimlmug. The restricted syndrome function erg depends only on the 
submatrix Hg. It is straightforward to see that the dimension of its image is the rank of 
Hg. □ 

Lemma 3.3. Let S be a stabilizer group of matrix H G M r ,n- The set Sg is an ¥2-vector 
space of dimension rank H — rank Hg- . 

Proof. The set Sg is the kernel of the F2-linear map: 

S — ► {E G V n I E C £} 

S I > S|£ 

By definition of the rank, the image of this application is a space of dimension rankHg, and 
the group S has dimension rankH. Therefore, dimSg = rankH — rankHg. □ 



From the lemmas, there are 2 rankH£ different syndromes and in each coset modulo S there 
are 2 rankH_rankH £ errors included in £. Therefore the number of correctable error patterns is 
2rankH+rankH £ -rankHg ^ an j s j nce there are 2 2 l £ l error vectors covered by £, the erasure vector 
£ can be corrected only if: 

2|£| < rankH + rank Hg - rank Hg. (4) 
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The rank of H is rank H = (1 — R)n where R is the rate of the quantum code. When p < 1 /2, 
there are typically more non-erased coordinates than erased ones, and it is reasonable to 
expect that the larger matrix Hg has a higher rank than the smaller matrix Hg. Equation ^ 
therefore becomes simply 2|£| < rankH and the typical weight of an erasure being |£| = np, 
we obtain: 

R < 1 - 2p 

which recovers ^ for the class of stabilizer codes. In the next section we will make this infor- 
mal argument rigorous and pave the way for improvements for particular classes of quantum 
codes. 



3.4 A combinatorial bound on the capacity 

Now, we will give a rigorous proof using an entropic formulation of this idea and Fano's 
inequality. 

Recall that a rate R £ [0, 1] is achievable if there exists a family of codes of rates Rt 
converging to R with vanishing error probability after decoding. Denote expectation by E. 
Let (Ht)t be a sequence of stabilizer matrices and denote by m the length of the code defined 
by the matrix Hj. 

Definition 3.4. The rank difference function D associated with the sequence {H-t)t of stabilizer 
matrices is D(p) = limsup t At(p) where : 

E p [rank H t g- rank H t £ ] 



n t 

Theorem 3.5. Achievable rates of a sequence of stabilizer codes of matrices (Ht)tgN> over the 
quantum erasure channel of probability p, satisfy: 

R<1- 2p- D{p). 

where D(p) is the rank difference function of the family (H.t)t- 

Proof. We shall apply the classical Fano inequality, see for example [15] . Recall that if X, Y 
and X are any three random variables such that X depends only on Y, then Fano's inequality 
states: 

log(|X|) 

where X takes its values in X. 

Over the quantum erasure channel, £ is the erasure vector random variable with distri- 
bution P(£ = v) = — p) n ~\ v \. The error random variable E is uniformly distributed 
among the errors acting on the erased components. We apply Fano's inequality when X is the 
information we need to recover the quantum state, namely the coset E.S of the Pauli error 
vector E. We set the variable Y to be the couple Y = (£, S) where £ is the erasure vector 
random variable defined above and S = cr(E) is the syndrome of E. The variable X is the 
best possible estimation of X given Y, meaning here that the decoding error probability P err 
is the probability to have 1^1. 

The conditional entropy is decomposed as: 

H(X\E, E) = £>((£, S) = (v, y))H(X\E = v, E = y). 

v,y 
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We have: 

H(X\E = v, E = y) = 2\v\ — rankH + rankH s — rankH,,, 

the proof of which is detailed in Lemma [3 . 6| below . We see that the value of H(X\8, = v, E = y) 
is independent of y, thus we have: 

H(X\E,Tt) = J^P(£ = t>)( 2 M - rankH- rank H„ + rank H c ) 

= 2np — rankH + E p (rankHg — rankHg). 

The random variable X = E.S takes on values in the quotient group X = V n /S. This 
quotient group is composed of |X| = 2 2n ~ rankH classes. From Fano's inequality we get, up- 
perbounding the denominator by 2n — rank H < 2n, 

2np — rank H + E p (rank Hg — rank Hg) — 1 

P< err — ^ 

In 

The rate of the quantum code is R = 1 — rankH/n. If the error probability goes to zero, then 
the rate of the quantum code family satisfies: 

limsupi? < 1 — 2p — D{p). 

□ 

Lemma 3.6. Let S be a stabilizer group of matrix H. The conditional entropy of X = E.S 
given £ = v and S = y is: 

H(X\E = v, E = y) = 2\v\ — rankH + rankH^ — rankH^, 

when the probability to have £ = v and E = y is non zero. 

Proof. Recall that given the erasure £, the distribution of E is uniform inside the support 
of £. Therefore, the probability of a coset E.S, assuming that the erasure is £ = v and the 
syndrome is E = y, is: 

WY FQ,c v ^ \{P€E.S\PCv,a(P) = y}\ 
F(X = E.S\E = v, E = y) = |{p - - { - - - <p) - ^ . 

When this probability is non-zero, by linearity, (multiply by an operator T C v of syndrome 
y), we can assume that y = 0. The set in the denominator is the subgroup N(S) V and the set 
in the numerator is a coset of the subgroup S v . Applying Lemmas 3J2 and |3,3| we have: 

I S" I 

PfX = E S\E = v E = y) = = 2 -2 l' u l +rankH-rankHs "'" rankHt ' 

\N(S) V \ 

whence 

H(X\E = v,Jl = y) = 2\v\ — rankH + rankH s — rankH,,. 

□ 



The following corollary proves the efficiency of our method. We recover the upper bound 
given by the capacity of the quantum erasure channel [5J- This bound is deduced only from 
combinatorial properties of stabilizer codes. It doesn't involve the no-cloning theorem. 
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Corollary 3.7. Achievable rates of a sequence of stabilizer codes of matrices (Ht)tgpj, over 
the quantum erasure channel of probability p, satisfy: 



R<l — 2p, 



when p < 1/2. 



Proof. To prove the corollary, it suffices to remark that A^(p) in Theorem 3.5 is non-negative 
when p < 1/2. Observe that we can write the function At as: 

A t (p) = cp t (i-p)-Mp), 

where 4>t{p) = ^(rankH^g)/?-^. It is intuitively clear and it is formally stated and proved in 
Appendix |A] (Proposition A. 3) that <pt is an increasing function of p. The corollary follows. □ 



Our goal is now to improve on Corollary 3.7 by finding non-zero lower bounds on D(p). 
This cannot be done for stabilizer codes in general since they are known to achieve capacity 
of the quantum erasure channel, but we can obtain such improvements for sparse quantum 
codes, i.e. codes that have sparse stabilizer matrices. Our most general result in this direction 



is Theorem 3.8 below. It is somewhat reminiscent of an upper bound on achievable rates of 
classical LDPC codes for the classical erasure channel |31) . 



3.5 Reduction to the study of the mean rank of a random submatrix of 
the stabilizer matrix 

Theorem 3.8. Let Q be any family of stabilizer codes of rates at least R and achieving van- 
ishing decoding error probability over the quantum erasure channel of erasure probability p. 
Suppose furthermore that every code C € C has a set of generators of its stabilizer group 
whose weights are all upper bounded by m. Then we have: 



R<(1- 2p) 



l-(l-p) 



m— 1 



l-(l-2p)(l-p) 



m— 1 ' 



Method: Let H be a stabilizer matrix and set <p{p) = E p (rankH£)/n. Denote by A the 
function A{p) = <p(\ — p) — (j){p). To apply Theorem 3.5 we need a lower bound on A(p). 



Because the function (f> is concave, any upper bound 4>(p) < M(p) implies the lower bound on 
A: 

. , , 1 — 2p ( rank H . , \ 

A(p) > P - M(p) . (5) 



1 — p \ n 

The formal proof of the concavity of (f> and of (j5j is somewhat technical and not necessary to 
the understanding of the main ideas, therefore it is placed in Appendix [A] (Proposition A.5). 



Proof of Theorem 3.8 Let H be a stabilizer matrix of a code C G C. Set 4>(p) = E p (rank Hg)/n. 



Let hg stand for the number of zero rows in the random submatrix Hg. We have: 



Hp) <~ [rankH-Ep(h^)] 

rankH rankH 

< (1-P) 



n 



n 
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Applying ^ we get: 



. , . 1 — 2p rank H . , _ 1 — 2p , „. . 

aCp > - — v - (i-p m >i — - i-p r 

1 — p n 1 — p 



From Theorem 13.51 we now have: 

R < 1 - 2p - (1 - R)(l - 2p)(l - p) m_1 

and the result follows after rearranging. □ 

As an example let us consider the family of color codes [7] . Color codes are defined from a 
trivalent tiling of a surface by faces and the associated stabilizer matrices have rows of weight 



bounded by the maximum length (number of edges) of a face. Hence Theorem 3.8 applies to 



this family that cannot be capacity-achieving if the faces stay with bounded length. 



4 Achievable rates of (2,m) CSS codes 

We now turn to deriving a refined upper bound on achievable rates of a particular family of 
quantum LDPC codes. Let us say that a binary matrix is of type (2, to) if every one of its 
rows is of weight m and every column is of weight 2. We shall say that a quantum code is a 
(2,m) CSS code if its stabilizer matrix H decomposes in two matrices Hx and H^, each of 
which is a (2, m) matrix. 

4.1 The 2-complex associated to a (2,m) CSS code 

The matrix Hj, viewed as a binary matrix, can be seen as the incidence matrix of a finite 
graph Gx- The vertex set V of the graph Gx is defined as the set of rows of Hy, and two 
vertices i and i! are declared to be incident if there is column j such that there are l's in 
positions (i,j) and The Edge set of the graph Gx can therefore be indexed by the 

columns of Hx The constant row weight to of means that the graph Gx is regular 
(every vertex has to neighbours). 

Recall that the classical code Cx is the set of vectors of F?; orthogonal to the rows of Hx- 
The code Cx is generated by the vectors whose supports coincide with cycles of the graph Gx 
(actually Cx is exactly the set of cycles of Gx when one allows cycles to be non-connected 
subgraphs) and Cx is classically called the cycle code of the associated graph Gx- 

If the graph Gx has n edges, the dimension of the code Cx is given by: 

Lemma 4.1. We have dimCx = n— \V\+kx where kx is the number of connected components 
ofG x . 

This is a classical result, see for example [6j. 

Since the rows of Hz are orthogonal to the rows of Hx, the supports of the rows of Hz 
are cycles of the graph Gx- The graph Gx together with the set of supports of the rows of 
Hz is called a 2-complex, and the supports of the rows of Hz are particular cycles that are 
called faces. That Hz is a (2, m) matrix means that every edge is incident to exactly two 
faces. 

In the same way that the matrix Hx defines a graph Gx, the matrix Hz defines a graph 
Gz, which together with Hx also makes up a 2-complex. The two complexes are said to be 
dual to each other: faces are vertices of the dual complex. 
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Figure 1: An example of self-dual 5-regular tiling of a surface of genus 4 composed of 16 
vertices, 40 edges and 16 faces. Each face is represented once. Vertices are represented several 
times to allow this planar representation. Each boundary edge is represented twice, multiple 
replicas of vertices and edges are identified to create the surface. In bold the identification of 
the edges {0, 1}. 

We shall say that the (2, m) CSS code (or equivalently the associated 2-complex) is proper 
if the two graphs Gx and Gz are connected and have girth (smallest cycle size) equal to m. 

It is not immediate that proper (2, m) CSS codes even exist, and the existence of families of 
(2, m) CSS codes with growing minimum distance is even less obvious. One way of coming up 
with such families is through the construction of combinatorial surfaces, the associated (2, m) 
CSS codes are a highly regular instance of surface codes [9|[37]. An example of such a surface is 
given on Figure[T]for m = 5. The associated (2, 5) matrices Hj and are given on Figure|2] 
It is the smallest (2, 5) CSS code we have found such that the associated graphs Gx and Gz 
are both connected and simple (without multiple edges). The only method we know of that 
allows the construction of proper (2, m) CSS codes involves sophisticated number-theoretic 
arguments and combinatorial surfaces. We shall take up this matter in Section [6] where upper 
bounds on the achievable rate of quantum (2, m) CSS codes lead to upper bounds on the 
critical probability for associated families of infinite tilings. 
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Figure 2: Two matrices of size 16 x 40 defining a (2,5) CSS code of parameters [[40, 10,4]]. 
Columns are indexed by the integers {0,2,..., 39} and rows are described by their supports. 
This code is the surface code defined from the 2-complex of Figure [T] 



4.2 Reduction to the study of the mean number of connected components 
in the subgraph of a graph 



Recall from our proof method described in Section 3.5 that our objective is to find an upper 
bound on the function (f)(p) = E p (rankH£)/n for a stabilizer matrix H. Since we are dealing 

with the stabilizer matrix H = |||^ ^ of a CSS code, we have <j)(p) = (f>x(p) + 4>z(p) where 

4>x(p) = IEp(rankHxe)/H and <f>z{p) = E p (rankH^£)/n. Recall that the two matrices Hx 
and are (2, m) matrices. Our problem is therefore to bound from above the mean rank of 
random submatrix of a fixed binary (2, m) matrix. 

In the rest of this section, let H £ M r ,n therefore stand for a binary matrix of type 
(2, m). We create the submatrix by keeping every column of H with probability p and 
independently of the others. 



The random subgraph G^: As described in section 4.1 we can regard the matrix H as 
the incidence matrix of a graph G with vertex set V. 

Given the random vector £ 6 Fj, we denote by C?£ the subgraph of G of incidence matrix 
i?£ . Assume that each component of £ is 1 with probability p and with probability 1 — p, 
independently. Then, the graph Gg is the random subgraph of G with unchanged vertex set, 
and obtained by taking each edge, independently, with probability p. In other words, taking 
a random submatrix of H corresponds to taking a random subgraph G^ of G. For the 



random subgraph G*£, Lemma 4.1 translates into: 



Lemma 4.2. If H is a binary (2, m) matrix with \V\ rows, then \V\ — ranking is equal to the 
number K£ of connected components of the subgraph Gg. That is: 

ranki?£ = |V| — Acg. 
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4.3 Bound on the mean number of connected components in the graph 



From Lemma 4.2 the mean rank of the submatrix can be expressed as a function of the 
expected number of connected components of the graph Gg: 

E p (rankff £ ) = \V\_ E p (^) 
n n n 

If we upper bound E p (ranking) /re by writing that E, p (k^) is lower bounded by the expected 
number of isolated vertices, we will simply recover Theorem 3.8 We will proceed to derive a 
more precise lower bound by enumerating larger connected components. 

Assuming that the m-regular graph G constructed from the matrix H has no small cycles, 
it looks like the m-regular tree G m in any sufficiently small neighbourhood. Thus, for our 
enumeration problem, we introduce the number of subtrees of G m , with k edges, containing 
a fixed vertex x of G m . We will use a generating function approach: for background, two 
classical references on this subject are |18] and |36| . 

The generating function for rooted trees: Let G m be the m-regular tree and let x 
be a fixed vertex of G m called a root. Let us define the generating function for rooted tree of 
degree m as the real function: 

T m {z) = ^2a k z k , 

fc>0 

where a k is the number of subtrees of G m , with k edges, containing the root x. To compute 
this generating function, it is useful to introduce the auxiliary generating function. 



k>0 

where bk is the number of subtrees T of G m such that : 

• T is composed of k edges, 

• the root x is included in T, 

• T contains a fixed edge {x, y} among edges incident to x, and no other edge incident to 
x is contained in T ■ 

This generating function doesn't depend on the particular choice of the edge x, y, by 
regularity of G m . The function T\{z) is sometimes called the generating function of planted 
rooted subtrees, since x has degree one in such a subtree. 

Coefficients bk can be computed easily using the Lagrange inversion Theorem |18] because 

satisfies: 

T 1 m {z) = z{l + T 1 m {z)r- 1 

This formula comes from the fact that every vertex of the tree except the root x has m — 1 
sons. We get: 

l/k(m-l) 
h= k{ k-1 

Then, the computation of the a k follows from the expression of T m as a function of 

T m (z) = (l + Z* (*))". 
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To see this formula remark that a subgraph of G m containing the root x can be decomposed 
into at m planted subtrees of root x. This method allows the computation of a large number 
of coefficients using symbolic computation software. 

We can now state an upper bound on the expected rank of the submatrix i?g involving 
the numbers a k - 

Proposition 4.3. If H is a binary (2, to) matrix whose associated graph G has girth (smallest 
cycle size) at least 5 + 2, then 

Ep(rankffg) 2 , ^m-2> 



< 

n to 



{i-(i- P ) m s 5 (p(i- P r- 2 )), 



where S$(z) = Y^k=o k^+l z ^ an< ^ ^ e a k are ^ e coefficients of the generating function T m . 

Proof. From ^ we want a lower bound on the expected number of connected components in 
the graph Gg. 

The graph Gg is constructed from the edge set of G by choosing each edge, independently, 
with probability p. Let us compute the expected number of connected components with 
edges. This is the average number of isolated points in the random subgraph Gg. Denote 
by Xq the random variable which associate with a random vector £, the number of isolated 
points in Gg. We can write Xq(£) = ^X„(£) where: X V (E) = 1 if the vertex v is isolated in 
Gg and X V (E) = otherwise. By linearity of expectation, we have: 

e(x ) =j2®( x v) = E p ( x - = x ) = iw -vT- 

V V 

Indeed, each vertex is bordered by to edges, therefore ¥(X V = 1) = (1— p) m , for every vertex v. 

This idea can be used with components of size (number of edges) k < S. For C a k-edge 
connected subgraph of G, let Xq denote the random variable equal to 1 if C is a connected 
component of the random graph Gg, and otherwise. The average number of connected 
components of size k is: 

E(x k ) = e ™ = vu^ 1 -p) m ^ 1 -p) m ~ 2 ) k - 

C connected 
subgraph of size k 



To prove the second equality we use two lemmas proved below. From Lemma 4.4, the ex- 
pected valu e of Xq is (1 — p) m (p(l — p) m ~ 2 ) k , independently of the subgraph G with k edges. 
Lemma 4.5 guarantees that the number of connected subgraph G is ^pjOfc. Finally, the quo- 

\v 



tient ^— is exactly This proves the proposition. □ 



Lemma 4.4. Let G be a m-regular graph of girth at least 5 + 2. If C is a connected subgraph of 
G with k < 5 edges, then C is a connected component of the random graph Gg with probability 
{l-p) m {p{l-p) m - 2 ) k . 

Proof. If k = 0, then G is an isolated point. This component appears in the random graph 
Gg with probability (1 — p) m by m-regularity of G. 

Assume the formula true for every connected subgraph C of size k — 1, with k < 5. A 
subgraph G of G of size A; — 1 is included in a ball of radius k — 1. We will prove that the 
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Figure 3: A ball of radius 3 in a 3-regular graph with girth > 7 

formula remains true if we add an edge to C. Let x be a vertex of C and let {x, y} be an edge 
which is not in C. Consider the graph C' = CL){x,y}. It contains k edges. Denote by dC the 
set of edges which have exactly one endpoint in C. Similarly, dC is the first neighbourhood 
of C . The set dC contains dC except {x,y}. Moreover it contains m—1 new edges: the 
edges {y, z} for z/i, These edges were not in dC. Indeed, if {y, z} is included in DC, define 
the graph C" = C U {y,z}. It contains k + 1 edges and it covers a cycle, because without 
the edge {y, z} it is still connected. This is impossible because the shortest cycle has length 
at least 5 + 2. Thus the formula is satisfied for all k < 5. □ 

Lemma 4.5. Let G be an m-regular graph of girth at least 5 + 2. The number of connected 
subgraphs of G with k < 5 edges is at least: 

|F| a 



jfe + 1 

where are the coefficient of the generating function T m . 

Proof. Connected subgraphs with k edges are trees and contain k + 1 vertices. It should be 
clear that given any vertex x, the number of ways of constructing a k-edge subgraph containing 
x, and hence the number of such subgraphs, is the same in G and in the m-regular tree. □ 

4.4 Achievable rates 

Let H = (hz) ^ e ^ ne stabilizer matrix of a (2, m) CSS code. Remark that the graphs 
associated to Hx and have girth at most m, since a row of Hz yields a cycle for Hx 
and vice versa. We will say that a (2, m) CSS code is proper if the two associated graphs are 
connected and have girth exactly m. 



We now translate the upper bound of Proposition 4.3 into a lower bound on the function 



D(p) of Theorem 3.5 for a sequence of (2, m) CSS codes. 

Proposition 4.6. For any sequence of proper (2,m) CSS codes we have: 

D(p) > (i^) (- - "(I - (1 -p) m Sm-2(p(l -P) m ~ 2 ))), 

VI— p J \m m 'J 

where S m -2{z) = ^2k=Q Fkl^ an< ^ ak are ^ e coefficients of the generating function T m . 
Proof. For H the stabilizer matrix of a proper CSS code, we apply ^ by using for the upper 



bound M(p) on E(rank(H£))/n, the sum of the upper bounds provided by Proposition 4.3 
on each of the terms E(rank(Hx,£))/fi and E(rank(H^g))/n. The result follows after some 
rearranging. □ 
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Finally, Proposition 4.6 together with Theorem |3. 5 1 yields: 



Theorem 4.7. Over the quantum erasure channel of erasure probability p, achievable rates of 
proper (2, m) CSS codes satisfy 



R < 



(1 - 2p)(— (1 - (1 -p) m 5 m _ 2 (p(l -p) m - 2 )) - l). 
\mp / 



where S m -2(z) = J2T=o Fkl 2 ^ an ^ a k are ^ e coefficients of the generating function T m . 



5 Erasure threshold of quantum LDPC codes 

In this part, we reformulate our bound on achievable rates of quantum LDPC codes by deter- 
mining an upper bound on the erasure decoding threshold of regular quantum LDPC codes. 

The capacity of the quantum erasure channel is 1 — 2p. This means that the rate of a 
family of quantum codes with vanishing decoding error probability over the quantum erasure 
channel of probability p, satisfies R < 1 — 2p. Alternatively, given a family of quantum codes 
of rate at least R, we can ask for the highest erasure rate that we can tolerate with vanishing 
error probability after decoding. This is the erasure decoding threshold. Assume that we 
have a family of quantum codes of rate higher than R, which achieve an asymptotic zero error 
probability over the quantum erasure channel of erasure probability p. Then, we have: 



P 



< 



1-R 



If we consider CSS codes of type (2,m), we have 2n/m rows in each matrix Hx and Hz- 
Therefore, the number of encoded qubits is at least (1 — 4/m)n (actually it is exactly (1 — 



4/m)n + 2/n by Lemma 4.1 ). Using the bound of Theorem 4.7 and the fact that the rates of 



the quantum codes are higher than 1 — 4/m, we obtain: 

Theorem 5.1. The erasure threshold of (2, m) proper CSS codes is below the solution of the 
equation: 



4 

m 



:i 



2 P ) ( — (i - (i - P rs m ^( P (i - _ i) 

\ mp v ' / 



where p G [0,1 /2] . 



This value is obtained at the intersection of the graphical representation of the upper 
bound with the line y = 1 — 4/m. An example of these curves is given in Figure [4] 

Using symbolic computation software, we computed different numerical values of this upper 
bound on the decoding erasure threshold Pe : 



type 


improved upper bound on p e 


Capacity bound p e < 2/m 


CSS(2,8), with Th. 


3.8 




0.228 


0.25 


CSS{2,8), with Th. 


5.1 




0.215 


0.25 


Stab(4,8), with Th. 


3.8 




0.228 


0.25 


CSS(2,5), with Th. 


3.8 




0.387 


0.40 


CSS(2,5), with Th. 


5.1 




0.381 


0.40 
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Th. 3.8 Th 4.7 l-2p 1-4/m 



rate R 
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0,1 0,2 0,3 0,4 

erasure probability p 



0,5 



Figure 4: The horizontal green line is the rate of a proper CSS (2,8) code. Its intersection with 
the capacity gives an upper bound on the erasure threshold of quantum codes : p e < 0.25. 



The two other curves are the upper bounds of Theorems 3.8 and 5.1 their intersection with 
the horizontal line gives the two upper bounds p e < 0.228 and p e < 0.215 for CSS (2,8) codes 
from Theorem 15. II 



6 Application to percolation theory 
6.1 Percolation theory 

In this section E denotes the edge set of a graph, rather than a Pauli error: context should not 
allow confusion. Let G = (V, E) be an infinite graph. Denote by [x v the probability measure 
on {0, 1} defined by /i p ({l}) = p. Consider the product space $7 = {0, 1} E endowed with 
the product probability measure P p = [i® E '. Random events should be seen as subgraphs. 
Informally, we choose every edge of G with probability p independently of the other edges, 
and obtain a random subgraph. The edges of this subgraph are called open edges. Percolation 
theory is interested in the probability that a given edge e is contained in a infinite open 
connected component (an open cluster). This probability depends a priori on the edge e, but 
not if the graph G is edge-transitive, for example if G is the infinite square lattice (Figure [5]). 
The central parameter in percolation theory is the critical probability p c , defined as: 



Pc (G) = inf{p G [0, 1], P p (|£(e)| = oo) > 0}, 

where £(e) denotes the open cluster containing edge e. 

By a famous result of Kesten |24j that stayed a conjecture for 20 years, we have p c = 1/2 
for the square lattice. Computing the critical probability exactly is usually quite difficult. 
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Figure 5: The square lattice 

For any integer m > 4, we denote by G{m) the planar graph which is regular of degree m 
and tiles the plane by elementary faces of length m. For m = 4 the graph G(4) is exactly the 
square lattice. The local structure of the graph G(5) is shown on Figure [6] 




Figure 6: The local structure of the graph G(5) 

For m > 4 these graphs make up regular tilings of the hyperbolic plane. Interest in perco- 
lation on hyperbolic tilings was raised in a number of papers e.g. [U [22] and determining 
their critical probability p c (rn) is highly non-trivial. Note that all graphs G(m) are self-dual 
like the square lattice G{m). 

We have the following easy bounds on p c : 

Proposition 6.1. The critical probability p c of G(m) satisfies 



Proof. We adapt the proof of [25] page 14 in the case of the square lattice. 

Let O be a fixed vertex. To show the first inequality we can say that there are not more 
than m(m — l) n_1 paths from O of length n in G(m) and the probability of such an open 
path is p n . So if p < j- the average length of an open path from O is not more than 
SnLi m ( m ~ l) n ~ 1 p n < oo. In this case p is under the critical probability. □ 
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The same method leads to the upper bound p c < 1 — The proof can be immediatly 

adapte from the case of the square lattice |25| page 14. 

6.2 Quotient graphs 

To study percolation on the hyperbolic tiling G(m), we need a family of increasingly big finite 
graphs which are locally the same as G(m). We will use a family introduced by Sirah in [33] . 

Let Pk(X) = 2 cos(fc arccos(X/2)) be the k-th normalized Chebychev polynomial and 
£ = 2 cos(-7r/m 2 ). Let y and z be the matrices of SX3 (Z[£]) defined by 

/ P m (0 2 ~ 1 P m (0 \ 

y = P m (0 1 

\ ~Pm(0 -1 / 




-1 -Pm(0 \ 

Pm{0 Pm{i?-l 
Pm(Z) Pm(0 2 1 / 



These two matrices generate the triangular group T(m) (3^ . To obtain a finite graph we can 
reduce the entries of the matrices modulo a prime number p. The coefficients are in the ring 
Z[£] which is isomorphic to the quotient X[X]/h(X) where h is the minimal polynomial of the 
algebraic integer £. Reducing coefficients modulo p, we obtain a group homomorphism from 
SX 3 (Z[f]) to SL 3 (¥ p [X]/(h(X)). The image of T(m) will be called f(m). 

Let G(m) be the graph defined like G(m) but with the group T(m), in other words the 
vertices, edges and faces of G(m) are defined as the left cosets of (y), (yz) and (z) respectively. 
There is a surjection s from G(m) to G(m) which sends u(y) to u(y). 

Following Sirah, let us define the injectivity radius of the graph G(m) as the largest integer 
r such that the restriction of the surjection s to a ball of radius r is one-to-one. It is shown in 
[33j that we can choose p so as to have r arbitrarily large. Loosely speaking, Sirah's argument 
is that if two distinct vertices u(y) and v(y) in G(m) have the same image under s then 
u~ 1 v in T{m) must project to the identity element in T[m). But this means that the matrix 
u~ 1 v has polynomial entries that, properly reduced modulo h(X), can only be expressed with 
coefficients at least one of which exceeds p: this implies that u~ 1 v can only be expressed as a 
product of a large number of matrices y and z, which in turn means that the original vertices 
u{y) and v(y) have to be far apart in G(m). 

The above construction enables us to define a family of finite graphs (Gy(m)) r >i such that 
each graph G r (m) has injectivity radius at least r, for every integer r. 

Let us now define random subgraphs of G r (m) through the product measure [J,p Er , where 
E r denotes the edge set of G r (m). In other words the open subgraph of G r (m) is created by 
declaring every edge open with independent probability p. 

For any fixed edge e, let £ r (e) be the (possibly empty) connected component of the random 
subgraph of G r (m) that contains e and call it again the open cluster containing e. Let f r (p) 
be the probability that |£ r ( e )| > r - We have: 

Proposition 6.2. If p < p c {m) then f r (p) goes to when r goes to infinity. 

Proof. Notice that the probability 1 — f r (p) that the open cluster containing e has cardinality 
not more than r is the same for the random subgraph defined on the finite graph G r (m) and 
the random subgraph defined on the infinite graph G(m). This is because this event depends 
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only on the ball of radius r centered on an endpoint of e, and these balls in G r (m) and G(m) 
are isomorphic. 

We can therefore consider f r (p) to mean the probability of the event F r that |£(e)| > r in 
the infinite graph G{m). Now (F r ) r >\ is a decreasing sequence of events, and ¥ p (D r >iF r ) is 
exactly the probability of percolation, which is since we have supposed p < p c . By monotone 
convergence we therefore have f r (p) = Pp(-F r ) —> 0. □ 

6.3 The quantum codes Q r (m) associated with the graphs G r (m) 

Every finite graph G r {m) gives rise to a CSS quantum code Q r (m) whose coordinate set is 
the edge set E of the graph. We will have therefore a quantum code of length n = \E\. The 
matrices Hx and Hz are defined as described in Section |4.1| The rows of are in one-to- 
one correspondence with the vertices of the graph. Every vertex x yields a row of Hx whose 
support is exactly the set of edges incident to x. Every row of Hx therefore has weight to. 
The rows of the other matrix Hz is in one-to-one correspondence with the set of faces of the 
graph. Every face yields a row whose support is equal to the set of edges making up the face. 
Since faces are m-gons, every row of Hz also has weight m. It should be clear that rows of 
Hx and Hz meet in either or 2 edges, so any row of Hx is orthogonal to any row of Hz 
and we have a quantum CSS code of type (2, to). 



Recall from Section 4.1 that the classical code Cx is the cycle code of the graph G r (m). 
When we reverse the roles of Hx and Hz by declaring the rows of Hz (rather than those of 
Hx) to be vertices and the rows of Hx to be faces, the graph thus defined is called the dual 
graph of G r (to) and we denote it here by G* (to) . The classical code Cz is thus the cycle code 
of the dual graph G*(m). 



From Lemma 4.1 since the graphs G r {m) and G*(m) are connected, we have: 



Proposition 6.3. The dimension k of the quantum code Q r (m) equals: 

k = ( 1- — J n + 2. 

\ TO J 

We remark that for m = 4, the graph G>(4) is a combinatorial torus and the quantum 
code Q r (4) is a version of Kitaev's toric code [26J. For to > 5 the quantum codes Q r (m) have 
positive rate bounded away from zero and minimum distance at least 2r (see the remark after 
the proof of Proposition |6.5| below) which is a quantity which behaves as logra. See |37] for a 
discussion of similar families of surface codes. 

Non-correctable erasures. An erasure vector can be identified with a set of edges of G r {m) 



(or of G*(m)) and we will denote it as before by £. From Proposition 3.1 we have that the 
erasure pattern £ is non-correctable if and only if £ either contains a cycle of G r {m) which is 
not a sum of faces (an element of Cx \ C%) or £, viewed as a set of edges of the dual graph 
G*(m), contains a cycle of G*(m) that is not a sum of faces of G*(m) (an element of Cz\C x ). 



6.4 Bounds on the critical probability using the capacity of the quantum 
erasure channel 

Consider an arbitrary member of the family of quantum codes Q r (m) associated with the 
graphs G r (m). Because the original graph G(m) is self-dual, all arguments involving G r (m) 
will be seen to hold for its dual graph G*(m) and we will focus on the probability that 
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the random erasure pattern £ contains a cycle that is not a sum of faces in the original 
graph G r (m). 

We would like to derive the upper bound on p c in Theorem |6.8| below by claiming the 
following: if p < p c , then for the family of graphs G r (m), the probability that the random set 



of edges £ contains a cycle which is not a sum of faces vanishes. If this is true, then Theorem 3.5 



applies and the rate R of the quantum code Q r (m) must satisfy R < 1 — 2p — D(p) for every 



p < p c and Proposition 6.3 gives the result since R = 1 — 4/m. 

Unfortunately, we do not know whether for every p < p c , the erasure pattern £ contains no 
cycle that is not a sum of faces with high probability. What we will prove however, is that if £ 
contains a cycle that is not a sum of faces, then with high probability one of the representatives 
of this cycle modulo the space of faces must have very small weight. To violate the capacity 
of the erasure channel we will therefore use, not Q r (m) directly, but an "improved" version 
Q' r (m) of Q r {m) that we now introduce. 

Proposition 6.4. Let Q r (m) be a hyperbolic code, n its length and R its rate. Suppose 
p e]0, ^[ and a G]0, 1[ are such that 

h(p) < q < ^, 

where h(p) = —p\og 2 p — (1 — /9)log 2 (l — p) denotes the binary entropy function. Then we 
can add an rows to the parity-check matrix Hx and an rows to the parity-check matrix Hz 
of Q r (m) to obtain a CSS code Q' r (m) of length n, rate R — 2a and distance d > pn. 

Proof. Denote by rx and rz the dimension of the code C x and C% respectively. We have 
rx = r z = - 1. 

We will construct a matrix H^- by adding an rows to the matrix Hx such that the rows 
of H' x are orthogonal to the rows of Hz an d the rank of is rx + cm. Let C' x be the code 
of parity-check matrix H^. 

For p G]0, l/2[, we define X p by 

X p (H' x ) = \{veC x \CMw(v)<pn}\. 
We can write X p as a sum a random variables to see that 

E(x) _ V l{H'> ^ C' x }\ 

v£B(0,pn) 

where B(0, pn) denotes the Hamming ball of radius pn. Let Li,L 2 , ■ ■ ■ L rx be rx rows of Hj. 
The number of suitable matrices is the number of families L'^L^,... L' an of vectors of ¥ 2 
such that V- G Cz for all j and (Li,L 2 , ■ ■ ■ L rx , L^L^ . . . , L' an ) are linearly independant. 

We can construct a suitable matrix if and only if rx + cm < dim(C^) this gives the 
condition a < (1 — — ) . In this case the number of matrices is 

rx+on-1 

JJ (2 n ~ rz -2 i ). 

i=r x 
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To evaluate the cardinality |{H^|t; € C' x }\ with v in Cx\C%, it suffices to add the 
condition U- S {y} for all j. We get 



rx+an—l 



So we have 



I! < 2 ' 

i=r x 

\{H! X \V G C' x }\ _ 2 n - r x- r z-*n _ -y 



< 2 



-a/1 



|{H^}| 2"-^"^ - 1 

This bound doesn't depend of v so we can give an upper bound on the expectation of X p 
because we know that the number of words in the ball of radius pn is less than 2 nh ^ . We 
find 

E(X P ) < 2 n ^- a \ 

If a > h(p) the mean goes to 0. Since X p has integer values there exists such that 
Xp(Hx) = 0. We obtain a CSS code of matrix with r' x = rx + an and Hz unchanged 
such that the minimum weight of a word of C' X \C Z is at least pn. 

We want to repeat this argument to have the minimum weight of a word of C' Z \C' X higher 
than pn. It suffices to choose a < ^(1 — + - because in this case rz + an < dim(Cx). □ 

Let £ be an erasure. We can write 

£ = £ c + £p (7) 

where £c is the sum of the connected components which do not cover a cycle which is not a 
sum of faces. The problematic part £p of £ is the union of the others components. 

In the graph G r (m), define g r (p) to be the probability that that the open cluster £ r (e) 
covers a cycle which is not a sum of faces. We have: 

Lemma 6.5. If p < p c (m) then g r {p) goes to when r goes to infinity. 

Proof. Recall that f r (p) denotes the probability that |£ r (e)| > r. We prove that g r (p) < f r (p) 



and apply Proposition 6.2 If |£ r (e)| < r then the open cluster £ r (e) is included in a ball of 
radius r of the graph G r {m). Since this ball is isomorphic to the ball of the same radius in 
the planar graph G(m), it is planar. In any planar graph every cycle is a sum of faces so £ r (e) 
covers a cycle which is not a sum of faces only if |£ r (e)| > r, hence g r {p) < f r (p)- D 

Remark: By the same planarity argument as above, every cycle of length less than 2r in the 
graph G r (m) is a sum of faces. This proves that the distance of the quantum code Q r (m) is 
at least 2r. 

Proposition 6.6. If we consider the erasure channel of probability p < p c then Ve > 0, 3ro 6 N 
such that if r > rQ then the expectation of the weight of £p defined as in ([7]) satisfies 

E(|£ P |) < en. 

Proof. For any edge e of G r (m), let X r ^ e be the random variable which take the value 1 if the 
connected component £ r (e) of e in G r (m) covers a cycle which is not a sum of faces and the 
value otherwise. Then we have: 

|£p| = ^2 X rie . 



To conclude note that E( J X" rje ) = g r {p) and apply Lemma 6.5 



□ 
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The next Lemma states that if the erasure vector £ has a large "problematic" part £ p then 



it must be correctable by the "improved" codes given by Proposition 6.4 



Lemma 6.7. Let Q' r (m) be one of the quantum codes given by Proposition \6.4\ and let d be 
its minimum distance. Suppose the part Ep of the erasure vector £ defined in ([7]) satisfies 
|£p| < d. Then £ is correctable by Q'(m). 

Proof. Denote by Cx and Cz the binary linear codes associated with the quantum code Q r (m) 
and by C' x , C' z their binary sub-codes associated with the quantum code Q' r {m) introduced in 



Proposition 6.4 and defined by augmenting the parity-check matrices Hx and Hz of Q r (m). 



If the erasure vector £ covers an element x of C^-\C^r then x must belong to Cx\C% 



i.e. 



x is a cycle of G r {m) which is not a sum of faces. The restriction of this cycle to £c defined in 
(J7| is another cycle y and the definition of Ec implies that y is a sum of faces. We obtain that 
x + y is included in Ep with y £ C% C C'^~, i.e. x + y S C' X \C'^- but this is a contradiction 
whenever the part Ep of the erasure £ has weight strictly less than the minimum distance d 
of the improved code Q'(m). □ 

We are now in a position to give an upper bound on the critical probability of G{m). Recall 



the definition of the rank difference function of a family of stabilizer codes (Definition 3.4). 



Theorem 6.8. Let m > 5, and let p c (m) be the critical probability for percolation on G(m). 
Let D(p) be the rank difference function of the sequence of stabilizer matrices associated to the 
tilings G r {m). Then for any p < p c we have: 

4 

1 < l-2p-D(p). 

m 



6.4 



Proof. Let R = 1 — ^ and fix p < p c . For any a such that < a < R/2, Proposition 
gives us a quantum code Q'{m) with minimum distance d> pn where p = h~ 1 (a/2) and rate 
R — 2a. For such a code the probability of a decoding error satisfies: 

Pew < P(\£p\ > pn). 



For any s > we can take r large enough so that Proposition 6.6 applies, and together 
with Markov's inequality we have 

Perr < P(\£-p\ > -En) < -. 

£ p 

For every e > we take p = ^/e. Then p(e) and simultaneously go to zero when e goes to 
zero. Defining a by a = 2h(p) and choosing a decreasing sequence of e's that tends to zero, 
we obtain a family of quantum codes Q' r (m) with decoding error probability tending to zero 
and rate R — 2a tending to R. 



We can therefore apply Theorem 3.5 to the sequence of stabilizer matrices of the codes 
Q' r (m). But we have just seen that their rates tend to 1 — 4/m and furthermore, since the codes 
Q' r (m) are obtained from the codes Q r (m) by adding a vanishing proportion of generators to 
their stabilizer group, the function D{p) is the same for the family Q' r (m) as for the family of 
surface codes Q r (m). □ 



Applying the lower bound on D(p) stemming from Proposition 4.6 we obtain after some 
rearranging: 
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Theorem 6.9. We have p c (jn) < p where p is the smallest solution p £ [0, 1] of the equation: 

i - ~ = a - 2 P ) ( — (i - (i - P rs m - 2 ( P (i - p) m - 2 )) - 1) . 

m V mp / 

where S m —2(x) = ^2™=$ Fkl 2 ^ an< ^ ak are ^ e coefficients of the generating function T m . 

Using symbolic computation software, we can compute this bound on the critical proba- 
bility. We use classical properties of generating functions to compute elements of the sequence 
(ofc)fe. For classical theorems on generating functions, in particular Lagrange inversion theo- 
rem, see for example |18| . |36j . 



m 


lower lower on p c (m) : j- 


bound of Th. El 


— = bound of Prop. 6.1 


5 


0.25 


0.38 


0.40 


10 


0.11 


0.16 


0.20 


20 


0.053 


0.073 


0.100 


30 


0.035 


0.046 


0.067 


40 


0.026 


0.033 


0.050 


50 


0.020 


0.026 


0.040 



We can observe that the new upper bound becomes better for tilings m with faces of large 
length. It is not surprising because, in this case, we enumerate the connected components of 
larger size. 

The exact value of the critical probability remains to discover. Numerical estimations of 
this value are difficult due to the exponential growth of balls of the graph. For example, Ziff 
pointed out the inconsistency of some numerical results in [22] • This reinforces the importance 
of our theoretical approach. 

7 Concluding Remarks 

• We have given a combinatorial proof of the upper bound on achievable rates of stabilizer 
codes R < 1 — 2p, over the quantum erasure channel. This proof is of course less general 
than previous proofs, since it applies only to stabilizer codes, but it gives mathematical 
insight into the degeneracy phenomenon and, as we have shown by the study of sparse 
stabilizer codes, it has the potential for improved upper bounds on achievable rates for 
particular classes of quantum codes. A generalization of this approach to the depolarizing 
channel would be very welcome. The difficulty is the fact that for depolarizing noise the 
probability of an error depends on its weight. We must deal with typical errors. 

• By graphical arguments, we proved that stabilizer and CSS codes defined by generators 
of bounded weight don't achieve the capacity of the quantum erasure channel. This 
result can be applied to surfaces codes, color codes and a lot of well studied families of 
quantum codes. This encourages us to construct irregular quantum LDPC codes and 
families with growing weights. 

• The exact value of the critical probability of hyperbolic tilings remains unknown. To im- 
prove our upper bound, we must enumerate more connected components in the random 
graph, which becomes more difficult when we count components that are not subtrees. 
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Our method also involves a concavity argument for the mean rank function which re- 
lates it to the rank difference function. This results in a manageable lowerbound on 
the relative rank difference function but it is generally not tight and alternative ways of 
evaluating this function would be desirable. 
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Appendix A: The rank of a random submatrix 

The goal of this section is to prove that the function 4>{p) = (rank He) is a concave (fl- 
convex) function. Then, we will use the concavity of cf> to obtain a lower bound on A(p) = 

</>(! -p) - 4>(p)- 

The key argument is the submodularity of the rank: 
Lemma A.l. Let H 6 Mr^iVi) be a Pauli matrix with n columns. The rank function: 

P({l,2,...,n}) — >N 

A i — > rank(vl) = rank(H^) 

is a submodular function. That is, the rank function satisfies: 

rank(A D B) + rank(yl Ufi)< rank(,4) + rank(B). 

This lemma embraces the case of a binary matrix. These rank properties can be obtained 
in the more general framework of matroid theory. A classical book on matroid theory is |29j . 

Proof. We will prove the lemma for binary matrices. Then, we will explain how to adapt our 
argumentation to the quaternary case. 

Let A and B be two subsets of {1,2, ... ,n}. From the dimension formula for the sum of 
two subspaces, we deduce: 

rank(^4 U B) = rank(A) + rank(5\,4) - dim(Imi^A) n (lmH B \ A ) 
rank(,4 U B) = rank(S) + rank(A\5) - dim(Imi^ s ) n (lmH A \ B ) 
rank(^4 n B) = rank(A) - rank(^4\5) + dim(Imi?AnB) n (lmH A \ B ) 
rank(^4 D B) = rank(S) - rank^^) + dim(ImiJ^ n s) n (lmH B \ A ) 

Look at the last term of these equalities. We have clearly ImH AnB C ImH A , therefore the 
space (lmH A f) B )Pi(ImH B \ A ) is a subspace of (lmH A )n(ImH B \ A ). This proves the dimension 
inequality: 

dim(lmH AnB ) n (lmH B \ A ) < dim(ImH A ) n (Im%). 

We have the same result exchanging A and B. We sum the four equalities and we apply the 
above inequality. We get the desired result: 

rank(A D B) + rank(yl Ufi)< rank(^) + rank(B). 
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To prove the property for a matrix with coefficients in V\, it is sufficient to show that 
all the tools from linear algebra used above are still satisfied for a Pauli matrix H. As 



recalled in Section 2.6 there is an isomorphism of F2 vector spaces between the Pauli group 
V n and the space F| n . Therefore, we can regard the matrix H 6 M T% rXP\) as a matrix 
[H X \H Z ] £ M r ,2n(^2)- The rank function can be written rankH = iank[H x \H z ] and the 



rank of a submatrix is: 

rankH £ = rank [H x \ iff ]. 
From this remark, the proof of the lemma is similar for stabilizer matrices. □ 



To study the derivatives of <f>, we introduce the function <3? depending on n variables 
x = (xi,X2, • • ■ , x n ) G [0, l] n defined by: 



®(x!,x 2 , ...,x n )= y~] 



(rankH £ ) J] 




This function can be seen as the expected rank E x (rank Hg ) for the probability measure such 
that the z-th component of £ is f with probability Xi and otherwise, independently of the 
other components. This polynomial function $ is infinitely derivable and its partial derivatives 
satisfy: 

Lemma A. 2. For all x in [0, l] n , we have: 

— (x)>0, Vie {1,2,..., n}, 

d 2 $> 

(x)<0, V*,ie{l,2,...,n}. 



dxidxj 

Proof. Let x be an element of [0, l] n . We remark that $ is an afline function in each variable. 
Thus, we have: 

— (x) = . . . , Xi-l, 1, X i+ i, ...,X n )~ $(xi, . . . , Xi-!, 0, x i+1 , ...,x n ) 

OXi 

= E rr (rank(£ U {*}) - E x (rank(£\{i}) 
= E X (^ rank(£ U {i}) - rank(£\{f} N 
> 0. 

In the preceding expression, the vector £ £ is considered as a subset of {1, 2, . . . , n}. Fixing 
Xi = 1 is equivalent to replacing the subset £ by £ U {i} and fixing x% = is equivalent to 
considering the subset £\{i}. 

Let j be an integer between 1 and n such that j 7^ i. The partial derivative §^:(x) is also 
an affine function of the j-th variable thus we can derivate it by the same process. We find: 

d 2 <$> , , d 9$, , 



dxjdxi dxj dx 



= E ;c (rank(£ U - E x (rank(£ U {i}\{j})) 

-E x (rank(£ U {j}\{t})) +E a .(rank(£\{» ) j})) 
= E x [rank(Aj U Aj) - rank(Aj) - rank(Aj) + rank(Aj n Aj)] 
< 0. 
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The subset Ai denotes the subset £ U and Aj denotes the subset £ U This 



□ 



quantity is negative by the submodularity of Lemma A.l 
If j = i then the second partial derivative is null. 

Proposition A. 3. Let H £ M rtn (Vi) be a Pauli matrix with n columns. The function (p(p) = 
^■E p (rankH£) is an increasing function on [0,1]. 

Proof. The function <f> is where i is the injection from [0, 1] to [0, l] n which sends p onto 

(p,p . . . ,p). The derivatives of <j> can be expressed as function of the partial derivatives of <£: 

^) = ^£p|!r(*(p))- 



1=1 



dxi 



From Lemma A. 2 this number is always non-negative. □ 

Proposition A. 4. Let H £ M rn (Vi) be a Pauli matrix with n columns. The function 4>(p) = 
^E p (rankH£) is a concave function on [0,1]. 



Proof. With the same notation as in the proof of Proposition A. 3 we obtain: 

d 2 § 



*"(p) = ~Ep 2 

n f—-* 



dxidxj 



(i(p)) 



□ 



From Lemma A. 2 this number is always non-positive, proving the concavity of 

Proposition A. 5. Let H e M rn (Vx) be a Pauli matrix with n columns. Assume that 4>(p) = 
^Ep(rankH) is upper bounded by M(p). Then, the function A(p) = 0(1 — p) — (j)(p) admits 
the lower bound: 



Hp) > 



1 - 2p\ /rankH 
1 — p J \ n 



M{p) 



Proof. It suffices to use the concavity of /. Indeed, by concavity the point (1 — p, 0(1 — p)) is 
above the segment between (p,4>(j>)) and (1,0(1)). That is: 



0(l-p)>0(p) + (l-2p) 



i — p 



The inequality follows using the equality 
M{p). 



rankH/n and then the upper bound (f>(p) < 

□ 
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